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We study the three-spin model and the Ising spin glass in a field using Migdal-KadanoiT approxi- 
mation. The flows of the couplings and fields indicate no phase transition, but they show even for the 
three-spin model a slow crossover to the asymptotic high-temperature behaviour for strong values 
of the couplings. We also evaluated a quantity that is a measure of the degree of non-self-averaging, 
and we found that it can become large for certain ranges of the parameters and the system sizes. 
For the spin glass in a field the maximum of non-self-averaging follows for given system size a line 
that resembles the de Almeida- Thouless line. We conclude that non-self-averaging found in Monte- 
Carlo simulations cannot be taken as evidence for the existence of a low-temperature phase with 
replica-symmetry breaking. Models similar to the three-spin model have been extensively discussed 
in order to provide a description of structural glasses. Their theory at mean-field level resembles the 
mode-coupling theory of real glasses. At that level the one-step replica symmetry approach breaking 
predicts two transitions, the first transition being dynamical and the second thermodynamical. Our 
results suggest that in real finite dimensional glasses there will be no genuine transitions at all, but 
that some features of mean-field theory could still provide some useful insights. 



I. INTRODUCTION 



Despite over two decades of work, the controversy con- 
cerning the nature of the ordered phase of short-range 
Ising spin glasses continues. Monte-Carlo simulations of 
three- and four-dimensional systems appear to be pro- 
viding evidence for replica symmetry breaking (RSB) in 
these systems (for recent reviews see 0J^)- However, 
recent developments have cast doubt on this interpreta- 
tion of the Monte Carlo data. In a series of papers on 
the Ising spin glass within the Migdal-Kadanoff approx- 
imation (MKA), we showed that the equilibrium Monte 
Carlo data in three and four dimensions that had been 
interpreted in the past as giving evidence for RSB can 
actually be interpreted quite easily within the droplet 
picture, with apparent RSB effects being attributed to a 
crossover between critical behaviour and the asymptotic 
droplet-like behaviour for small system sizes We 
also showed that system sizes well beyond the reach of 
current simulations would probably be required in order 
to unambiguously see droplet-like behaviour. Very re- 
cently, a third view on the nature of the low-temperature 



phase if spin glasses has emerged. In this picture, there 
exist droplet excitations in short scales, but on large 
scales there are system-wide excitations that cost only a 
finite energy in the thermodynamic limit, and that have 
a surface whose fractal dimension is less than the space 
dimension |p|-p^. It remains to be seen whether these 
excitations survive at larger system sizes. Within MKA, 
fractal excitations are not possible, and the signatures of 
these excitations found in Monte-Carlo simulations are 
therefore not present in MKA. 

There is a close connection between the question of the 
nature of the spin glass phase and that of the existence 
of a phase transition in a spin glass in an external field. 
Mean-field theory predicts a phase transition to a spin 
glass phase with RSB along the so-called de Almeida- 
Thouless line. The droplet picture predicts no transi- 
tion. The reason is that in the presence of a field time 
reversal symmetry is broken, and there is no symmetry 
left that could possibly be broken by a phase transition, 
except for replica symmetry. Monte-Carlo simulation of 
the spin glass in a field ||l^-|l5|] show some evidence of a 
phase transition, and in particular of non-self-averaging, 
however, the situation is complicated by the presence of 
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large finite-size effects due to crossover phenomena. For 
this reason, Parisi et al jl6| studied a different system 
that also has broken time-reversal symmetry, but is ex- 
pected not to have strong crossover effects. This system is 
the three-spin model, where the two-spin products in the 
Ising spin glass without field are replaced by three-spin 
products. The numerical evidence for a phase transition 
in the four-dimensional system seems good, if a quantity 
that measures the degree of non-self-averaging is studied. 

It is the purpose of this paper to study the two men- 
tioned systems without time-reversal symmetry in MKA, 
in order to check whether a similar degree of non-self- 
averaging could be produced by finite-size effects. The 
flow of the coupling constants shows that for both sys- 
tems there exists only one attractive fixed point which 
corresponds to a paramagnet in a random field, and that 
there is consequently no phase transition within MKA. 
Nevertheless, in both systems the coupling constants in- 
crease initially for sufficiently low temperatures, indicat- 
ing that for small system sizes there might be the appear- 
ance of a phase transition. Then, we looked at the non- 
self-averaging parameter in both models within MKA, 
for various system sizes and parameter values. We found 
a behaviour similar to the one reported for the Monte- 
Carlo simulations, and apparent RSB for system sizes 
similar to theirs. Furthermore, for the spin glass in a 
field, the maximum of the non-self-averaging parameter 
as function of the field (for fixed system size) marks a line 
that can be interpreted as a remnant of the de Almeida- 
Thouless line. 

Some insights into what might be expected in the fi- 
nite dimensional three-spin model can be obtained from 
the mean-field solution of the p-spin model (for a re- 
view of which see [0). It has an analytical solution in 
the limit where p goes to infinity which can be obtained 
by a one-step replica symmetry breaking scheme. For 
our purposes the solution is best understood in terms of 
metastable states which can be identified with the solu- 
tions of Thouless- Anderson-Palmer (TAP)-like equations 
psf . The partition function is obtained from the integral 

Z= [^"^^\feMNa{f)-NfP), (1) 

where P = l/ksT and N is the number of spins in the 
system. / is the free energy density of a TAP state. So- 
lutions of the TAP equations exist for fmin < f < fmax, 
and the number of solutions at free energy density / is 
exponentially large and given by exp(A^cr(/)), with (t(/) 
vanishing at both fmin and fmax- For temperatures in 
the interval < T < T5 the integral is dominated as 
^ 00 by the lower limit of the integral, i.e. states 
whose free energy Nf differs from the state of lowest en- 
ergy by only a finite amount. For temperatures in the 
interval T5 < T < the integral can again be done by 
steepest descents and is dominated by some value of / ly- 
ing in the interval fmin < f < fmax- As the temperature 
approaches To this value tends to fmax- In the tem- 



perature interval Ts < T < Tjj an exponentially large 
number of states contribute to the thermodynamics in 
contrast to the situation below Ts where only a finite 
number contribute. Above T^, only the trivial para- 
magnetic state contributes. There are thus two phase 
transitions at mean-field level. The lower temperature 
transition at Ts is accompanied by singularities in the 
free energy but at the higher temperature transition T^ 
the free energy is smooth and the presence of a transition 
is best inferred from singularities in the dynamics. 

Now for a finite dimensional system metastable states 
with / > fmin are unstable. (Imagine in such a state 
converting a block of spins of linear dimension L to have 
the orientations that they would have in the lowest state; 
the free energy gain will be or order L''; the energy cost 
of creating such a region will be no more than the cost 
of breaking all the bonds at the surface of the region, 
^d-i_ fpjj^^g possibility of nucleating lower free energy 
states prevents the existence of metastable states in finite 
dimensional systems). As the transition at the higher 
temperature To involves the metastable states only, one 
deduces that it will not exist in a finite dimensional sys- 
tem. The only transition which could possibly survive 
to finite dimensions is the one associated with Ts- Our 
studies however indicate that it too probably does not 
occur in finite dimensional systems. 

Our study of the three-spin model had another, per- 
haps physically more significant motivation. The three- 
spin model and its cousins have been extensively stud- 
ied at mean-field level as models of structural glasses 
]l9[ . The higher temperature transition, Td, whose sig- 
nature is purely dynamical (and very similar to that in 
the mode-coupling approach to real glasses [^), while 
the transition at the lower temperature, Ts, is associ- 
ated the Kauzmann temperature, Tk the tempera- 
ture at which the configurational entropy of the glass goes 
to zero. It has been a common belief of many workers 
for several decades that there is no genuine transition at 
Tk- Recently this belief has been strongly reinforced by 
the novel Monte Carlo simulation of Santen and Krauth 
|p2| , who found no evidence of a genuine transition in a 
simulation in which the particles could be properly equi- 
libriated. Our work strengthens the argument that no 
genuine transition analogous to Ts or Tk will exist in 
finite dimensions. Nevertheless we can see echoes of the 
mean-field results in our calculations. It is our belief 
that future work should focus on how the singularities in 
mean- field results are " rounded-ofF' in finite dimensions. 

This paper is organized as follows: First, we define the 
models and quantities studied in this paper. Then, we 
describe the MKA for the two models studied in this pa- 
per. We use two different methods in order to make sure 
that the results do not depend on the particular imple- 
mentation of the MKA. Next, we study the three-spin 
model and give our results for the flow of the coupling 
constants, and for the non-self-averaging parameter. In 
section V, we discuss the spin glass in an external field. 
Again, we give results for the flows and the degree of non- 
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self-averaging. We also give scaling arguments based on 
the droplet picture that explain most of our findings. Fi- 
nally, we summarize and discuss our results. 



II. MODELS AND DEFINITIONS 

The Edwards-Anderson spin-glass Hamiltonian H in 
the presence of a uniform external magnetic field h is 
given by 



(2) 



where the Ising spins <7i can take the values ±1, and 
the nearest-neighbor couplings Jij are independent from 
each other and Gaussian distributed with a standard de- 
viation J. For convenience, the Boltzmann-factor (3 is 
absorbed into the couplings and fields. Without a field h, 
the model has a low-temperature phase with nonvanish- 
ing correlations (SiSj) even for spins that are far apart. 
According to the droplet picture, this phase is unique (up 
to a global flip of all spins), and it is destroyed as soon as 
the field is turned on. The reason is that the field induces 
regions of a sufficiently large radius to flip if the magne- 
tization of this region opposes the field. The radius r 
of these regions is obtained from the condition that the 
gain in magnetic energy, hr'^^'^ becomes comparable to 
the loss in coupHng energy, Jr^, leading to 



(J//i)i/('^/2-«). 



(3) 



Here, d is the dimension of the system, and 9 the scal- 
ing dimension of domain walls. Beyond the radius r, 
the long-range correlations of the spin-glass phase are 
destroyed. 

In contrast, the RSB picture predicts the existence of 
infinitely many different phases of comparable free en- 
ergy in the absence of a field. With increasing field /i, 
the number of phases decreases, and it becomes one at 
the de Almeida-Thouless line hc{J). At the critical spin- 
glass transition J^, the critical field vanishes, and it 
diverges to infinity as the coupling strength J diverges 
(i.e., as the temperature goes to zero). 

The numerical analysis of the spin glass in a field is 
hampered by strong crossover effects. Crossover effects 
are expected to be less strong in the three-spin model, be- 
cause it has no tunable parameter that can restore time 
reversal symmetry and lead to strong crossover effects 
when small [|l5|. In the most easily tractable version of 
this model, each site is occupied by two Ising spins, ct|^'' 

(2) 

and , and the Hamiltonian is given by 



+ ^VWa«af +J^)af)a«.f), (4) 



where ij are nearest-neighbor pairs, and the couplings 
Jj^"^ are chosen independently from a Gaussian distribu- 
tion with zero mean and width J. When the signs of all 
spins are reversed, the sign of the Hamiltonian changes 
also, indicating the violation of time-reversal symmetry. 

If finite-dimensional systems have no RSB, this model 
has no phase transition since there is no symmetry that 
could be broken. On the other hand, if RSB occurs 
in finite-dimensional spin glasses, the three-spin model 
could show a phase transition at some critical coupling 
strength Jc. 

It has proven useful to consider two identical copies 
(replicas) of the system, with the spins {(Ti} and {t^}, 
and to measure overlaps between them. This gives infor- 
mation about the structure of the low-temperature phase, 
in particular about the number of pure states. The main 
quantity studied in this paper is the parameter A which 
measures the degree of non-self averaging, and is defined 

by 



A 



[(('Z-('Z»= 



2\21 



(5) 



where (...) and [...] denote the thermodynamic and disor- 
der average respectively. The overlap q between the two 
replicas is given by 



for the Ising spin glass in a field, and by 



q = 



for the three-spin model. N is the number of sites in the 
system. 

A is most easily evaluated by introducing a coupling 
between the two replicas, and by differentiating with re- 
spect to it. The Hamiltonian for the coupled system is 
then 

-/3K(a, r) = -/37^((7) - /37^(t) eiVq. 
The mean overlap is given by the expression 



and its variance by 

(('Z-('Z»') 



l^lnz' 

N dt 



1^ 

TV a? 



InZ 



e=0 



Z is the partition function. 

In systems with RSB, the probability distribution P((i) 
of q is broad, and A has a nonzero limit in the limit of 
infinite system size. On the other hand, in the absence 
of RSB, each sample has only one, sample-independent 
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value of gr, and A vanishes in the thermodynamic hmit. 
Consequently, an A that increases with increasing system 
size could be taken as an indicator of RSB. However, we 
will see in this paper that even systems without RSB can 
show an increasing A over a wide range of system sizes 
and parameters. 



III. THE MIGDAL KADANOFF 
APPROXIMATION 

The Migdal-Kadanoff approximation (MKA) is a real- 
space renormalization group that gives approximate re- 
cursion relations for the various coupling constants. An 
exact decimation, which consists in taking the trace over 
all those spins that do not belong to the coarse-grained 
lattice, generates higher-order couplings between spins of 
more that two sites, and is therefore untractable. In order 
to circumvent this problem, the MKA moves the bonds of 
a hypercubic lattice before each decimation step in such 
a way that no higher-order couplings can be generated. If 
the bond-moving shall be symmetric with respect to the 
different space directions, one ends up with the scheme 
represented in Figure |l|. 




FIG. 1. The Migdal-Kadanoff bond-moving scheme for a 
cubic lattice. 



In a d-dimensional lattice, 2'^^^ bonds are superim- 
posed as a consequence of bond-moving. In the absence 
of field terms (i.e. of terms that couple only spins that 
sit on the same site «, and are therefore not in a clear 
way associated with bonds), the 2'^~^ coupling constants 
of each of the d bundles of bonds per coarse-grained unit 
cell simply add up, and the "naked" spins that are left be- 
hind have no couplings. Taking the trace over the d spins 
that sit on the d main bonds leads to coarse-grained cou- 
pling constants between neighboring spins on the coarse- 
grained lattice. Taking the trace over the naked spins, 
gives only a constant contribution to the partition func- 
tion, which can be neglected. This decimation procedure 
is iterated n times on a lattice of linear size L = 2", until 
a single unit cell is left over. Assuming periodic bound- 
ary conditions, one can then take the trace over the final 
spin. 

The flow of the coupling constants in this scheme re- 
sults from alternating the addition of 2'^~^ bonds with 
linking two of these new bonds together and taking the 
trace over the middle spin. Essentially the same flow re- 
sults when a decimation is done on a hierarchical lattice 
that is constructed iteratively by replacing each bond by 



2^^ bonds, as indicated in Fig. ||. The total number of 
bonds after n iterations is 2'*". Spin decimation on such 
a lattice is done by taking the trace over the spins that 
are highest on this hierarchy (i.e., that were added last 
during the construction procedure). At each decimation 
step, first the trace is taken over the middle spin of two 
linked bonds, and then 2'^~^ bonds are added together 
to form a new bond, until the lowest level is reached 
and the trace over the remaining two spins is calculated 
. Apart from the fact that the order of bond-adding 
and decimation is reversed, the recursion relations for the 
coupling constants are obtained by the same procedure 
as for the bond-moving algorithm. 




FIG. 2. Construction of a hierarchical lattice. 

This equivalence between a bond-moving procedure for 
a hypercubic lattice and a hierarchical lattice hold no 
longer when field terms are present. In both of our mod- 
els, these field terms are present either from the begin- 
ning (Ising spin glass in a field), or they are initially ab- 
sent, but are generated during the decimation procedure 
(three-spin model). For the bond-moving procedure, one 
has to decide whether the fields shall also be moved (and 
to where), or whether they shall remain at their sites, or 
whether part of them shall be moved. This creates a cer- 
tain freedom in the renormalization scheme, and the most 
plausible choice is determined by the requirement that 
the fiows of the fields near the zero-temperature fixed 
point and at the infinite-temperature fixed point shall 
be those of the hypercubic lattice. However, if the main 
results of the MKA shall be generic, they should not de- 
pend on the precise implementation of the bond-moving 
algorithm. Otherwise one might doubt that the MKA 
reflects the features of the real system. For this reason, 
we have performed the MKA for a variety of different im- 
plementations, and we found our main conclusions con- 
cering A to be independent of the implementation. In 
the following sections, we will give our results for the hi- 
erarchical lattice and for a bond-moving scheme that has 
the correct flow of the field terms. 

The treatment of a hierarchical lattice in the presence 
of field terms is straightforward. In order to understand 
that the flows of the fields are different on a hierarchical 
lattice compared to a hypercubic lattice, let us consider 
a situation where the fiows of the couplings go to zero 
with increasing iteration number, which is the situation 
that we will encounter below for both models. As long 
as the couplings are nonzero, each decimation step gen- 
erates a contribution to the fields at the sites that are 
left over. The two corner spins, which are left over until 
the end, receive consequently 2"'^'^^^^ field contributions 
from the first decimation, 2^"^^^''^^^^ from the second it- 
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eration, and so on, until the couplings are virtually zero. 
For twice the system size, i.e., for a lattice with n + 1 lev- 
els, the mean of the field contribution to the corner spins 
due to decimations is larger by a factor 2^'*"^^ and so is 
the variance of the field contribution. Even though the 
couplings go to zero after a certain number of iterations, 
the fields keep growing. 

In contrast, on a hypercubic lattice, the fields must re- 
main constant as soon as the couplings have become zero. 
Clearly, this can only be achieved if fields terms are not 
moved to the sites that will not be traced over. On the 
other hand, field terms must not be left with the "naked" 
spins. The reason is that near the zero-temperature fixed 
point where the couplings are very large, all fields must 
add up under renormalization. Fields must therefore al- 
ways stay with spins that are coupled to other spins. For 
this reason, fields should be moved to those d sites that sit 
at the middle of the main bonds. Even with this restric- 
tion, there remains some freedom in choosing which field 
should be moved where. In our simulations, we treated 
field terms as belonging to bonds. The initial fields were 
evenly distributed between the ends of all bonds, and the 
fields generated during decimations naturally end up at 
the ends of bonds. When a bond was moved, we moved 
all its field terms to that end that was to be traced over 
next. 

For the Ising spin glass in a field, the recursion of four 
different parameters must be considered when studying 
the flow diagram and thermodynamic quantities. These 
are the two-spin coupling, the two fields on the two ends 
of a bond, and a constant. If one evaluates quantities 
related to the overlap between two replicas, each site has 
two spins, leading to 16 parameters. The same number 
of parameters occurs for the three-spin model, if only one 
replica is needed, as, e.g., for the flows and the phase di- 
agram. For the evaluation of A, we need two replicas, 
leading 256 couplings. Luckily, the decimation step can 
treat all 256 parameters with the same formula, which 
involves a 256 x 256 matrix that is calculated once at the 
beginning of the program. 



IV. THE THREE-SPIN MODEL 

This model was studied using Monte-Carlo simulations 
in four dimensions in Jl5| , ^,0^ , and evidence for RSB 
was found. The authors of P3[p^| found in particular that 
the non-self-averaging parameter A is small for large tem- 
peratures, and becomes large for smaller temperatures. 
Curves for different system size L = 3, 4, 5, 6 cross nearly 
at the same temperature, and below this temperature A 
increases with increasing L. Thus, the degree of non- 
self averaging increases with the system size, just as can 
be expected for a replica-symmetry breaking transition. 
Monte-Carlo simulations ||lj,^,|l6| arc usually done with 
couplings J — ±1. The precise distribution of the cou- 
plings should however not affect the universality class. 



Analytical results were obtained for the p-spin model in 
mean-field theory, where one-step RSB was found. This 
means that the ground state has a nonzero probability of 
being occupied below a critical temperature Ts (see Sec- 
tion W . This mean-field scenario is fundamentally differ- 
ent from the full RSB claimed to be seen in Monte-Carlo 
simulations of the four-dimensional system. Thus, the 
argument usually employed for spin glasses that mean- 
field like behaviour can be found in finite-dimensional 
short-range systems, fails here. 

In the following, we show using MKA that the assump- 
tion of the absence of crossover effects in this model is 
incorrect, and that A might at low temperatures and for 
small system sizes increase with increasing system size 
even if the system is self-averaging in the thermodynamic 
limit. We mainly focus on the case of four dimensions, 
but report also some results in d=2 and 3. Let us first 
discuss the flow of the coupling constants as the system is 
renormalized. Because each bond is connected to 4 spins, 
the flow of 16 coupling constants has to be considered. In 
order to obtain this flow, we iterated the recursion rela- 
tion on a set of 10^ bonds. At each iteration, each of the 
new set of 10® bonds was generated by randomly choos- 
ing 2^ bonds from the old set. For a hierarchical lattice, 
where the generated fields remain at that end of a bond 
at which they are generated, we first took the trace over 
the inner spins of each of the 2'^"-^ pairs of bonds, and 
than we added the resulting bonds; for the bond-moving 
procedure described in the previous section, we first gen- 
erated two bunches of 2"^^^ bonds each, then moved the 
fields of all but the "original" bond of each bunch to the 
inner spin, and took the trace over the inner spin. Figure 
^ shows the flow of the width of the three-spin couplings 
for different initial width in four dimensions, for the two 
different algorithms. 




2 4 6 8 



iteration 

FIG. 3. Flow of the three-spin coupling strength J for 
bond-moving (top) and the hierarchical lattice (bottom) for 
d — 4, divided by the initial coupling strength. The 
curves correspond to J = 0.065, 0.13, 0.18, 0.26, 0.5552, 1, 2, 5 
(bond-moving) and J = 0.13,0.26,0.35,0.52,1,2,5,10 (hier- 
archical lattice) from bottom to top. 
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One can see that for weak coupling (or, equiva- 
lently, high temperature) the coupHng strength de- 
creases quickly with increasing system size. However, 
for stronger coupling or lower temperatures, the coupling 
strength increases during the first few iterations, and de- 
creases afterwards. The maximum is reached between the 
3rd and 4th iteration, or between L = 8 and L — 16. For 
sufficiently strong coupling, the curves reach an asymp- 
totic shape. On the hierarchical lattice, where the fields 
grow without bounds, the 3-spin couplings decrease to 
zero faster than with bond-moving. Furthermore, curves 
for the hierarchical lattice seem to correspond roughly 
to those of the bond-moving procedure if the three-spin 
couplings are divided by a number around 3. The reason 
is that the first step during the bond-moving procedure 
summarizes 8 bonds to one new bond. The width of the 
three-spin coupling is therefore increased by a factor of 
\/8 in four dimensions. In order to compare to the hi- 
erarchical lattice or to Monte-Carlo simulations on a hy- 
percubic lattice, one should divide the coupling strength 
of the bond-moving procedure by VS- 

If one considered only systems of sizes up to 8, one 
would get the illusion of a phase transition with a (1/ J)c 
around 3 or 4, a value which is not far from the one given 
for Tc in |jl^ . (Note that these authors have kept the cou- 
pling strength fixed at ±1, and varied the temperature. 
Their T corresponds therefore to our 1/J.) 

Figure ^ shows the flow of the widths of the different 
coupling constants for an initial width of the three-spin 
coupling J = 2. 




FIG. 4. Flow of the width of the field (long dashed), of 
the on-site two-spin coupling (dotted), the two-spin coupling 
across a link (dashed), the three-spin coupling (solid), and 
the four-spin coupling (dot-dashed) for bondmoving (top) and 
on the hierarchical lattice (bottom) for d = 4, for an initial 
three-spin coupling J = 2. 

One can see that the strengths of the field and of the 
on-site two-spin coupling (which can also be viewed as a 
"field" ) increase rapidly and without limits for the hier- 
archical lattice, and that they saturate at a finite value in 
the bond-moving case. The two- three- and four-spin cou- 



plings increase during the first few iterations, and then 
decrease again. Thus, our three-spin model corresponds 
on large scales to a system only with random fields and 
random couplings between the a and t spins on the same 
site. There are no couplings between spins on different 
sites on large enough scales, but sites are independent 
from each other. Only on small scales could one get the 
impression that the system has long-range correlations. 
However, these system sizes are exactly the ones studied 
in |l5|l22 



The crossover regime becomes larger with increasing 
dimension. Figure ^ shows the flow of the three-spin 
coupling for an initial value J = 10 in d = 2, 3, 4 dimen- 
sions. Clearly, the strength of the increase and the range 
of system sizes over which this increase occurs increases 
with increasing dimension. One can therefore expect that 
in even higher dimensions, the apparent phase transition 
becomes more pronounced. 




FIG. 5. Flow of the width of the field for an initial value 
J = 10 in d = 2,3,4 dimensions (from bottom to top). The 
top graph is for bondmoving, the bottom graph for a hierar- 
chical lattice. 

Next, let us study the non-self-averaging parameter A. 
As explained in the previous section, A can be evaluated 
by introducing a coupling between two identical replicas 
of the system. Since there are now 8 spins associated 
with each bond, the number of couplings that have to 
be evaluated in MKA approximation is 2^ = 256. Fig- 
ure ^ shows A as function of the coupling strength for 
different system sizes up to 16, in 4 dimensions. Larger 
system sizes could not be studied due to limitations in 
computer time. One can see that A increases with in- 
creasing system size whereever it is appreciably different 
from zero, and reaches large values. This figure gives 
the impression that the system shows non self-averaging. 
Of course, for larger system sizes, A must eventually de- 
crease again since we know that there is self-averaging 
in the thermodynamic limit. In contrast to the Monte- 
Carlo simulation results |23 1^, our curves for A do not 
intersect at a coupling strength and A value of the order 
1. We have performed a similar simulation in d = 2 di- 
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mensions and found that A increases as the system size 
increases over the range L = 2,4,8,16. However, for 
L = 32 and L — 64, A decreases. If we assume that 
the system size for which A is largest increases with each 
dimension by a factor 2, as it does for the flow of the 
couplings, we can expect that in rf = 4 the system sizes 
for which a decrease in A can be seen is beyond L = 64. 




FIG. 6. The non-self-averaging parameter A 
for L = 2,4,8,16 (from bottom to top), and d = 4. The 
average is taken over 10000 samples for the smallest system 
size and 100 samples for the largest. The top graph is again 
for bond-moving, a nd the bottom graph on the hierarchical 
lattice. 

To summarize this section, we have shown that the 
three-spin model, even in situations where we know that 
it self-averages in the thermodynamic limit, can show 
indications of non-self-averaging at those system sizes 
typically studied in simulations. Evidence for non-self- 
averaging found in Monte-Carlo simulations must there- 
fore be taken with caution as it might be misleading. 



V. ISING SPIN GLASS IN A MAGNETIC FIELD 

Monte-Carlo simulations in four dimensions show some 
indication of RSB [p"3|-p^. Just as for the three-spin 
model and for the spin glass without external field, these 
findings may again be due to finite-size effects and to the 
closeness to the critical temperature. Indeed, an inves- 
tigation of the ground-state structure of a spin glass in 
a magnetic field |2^ shows no indication of RSB. (See, 
however, the discussion in p^ , p7t .) 

In order to test for finite-size effects, we studied the 
spin glass in a field using MKA, and determined the 
non-self-averaging parameter A as function of the sys- 
tem parameters. We found that the degree of non-self- 
averaging can be large for the system sizes typically used 
in simulations, in particular when the contribution of the 
field to the free energy is comparable to that of the cou- 
plings. While most published Monte-Carlo simulations 
were done in four dimensions, we chose to study the MKA 



in three dimensions, in order to be able to go to larger 
system sizes. Because there are three parameters to be 
varied (the system size, the field, and the two-spin cou- 
plings), many data points had to be collected, and this 
is done faster in 3 dimensions. Of course, we expect that 
the results of the MKA are similar in four dimensions, if 
the exponents for 3 dimensions are replaced with those 
for 4 dimensions. Just as for the three-spin model, the 
apparent non-self-averaging should become even stronger 
in 4 dimensions. 

First, let us study the flows of the couplings and fields. 
The decimation procedure leads to the creation of ran- 
dom fields, while the mean value of the field is not 
changed. Figure |^ shows the flow of the two-spin coupling 
J for various initial values, and for a flxed field h — 0.1. 
For initial couplings larger than the critical coupling (1.13 
for the hierarchical lattice and 0.55 for bondmoving), the 
coupling strength decreases immediately. However, if the 
initial coupling strength is sufficiently deep in the low- 
temperature phase, it increases first, until the random 
field has become strong enough to have a reducing effect 
on the coupling strength. Ultimately, the fiow goes to a 
fixed point where the coupling strength is zero. On the 
hierarchical lattice, the width of the field keeps growing 
indefinitely, while it saturates in the bond-moving case, 
as discussed in section lU. Clearly, there is no phase 



transition in the presence of an external field, but there 
are strong crossover effects if the field is small. 




FIG. 7. Flow of the two-spin coupling strength J, divided 
by the initial strenggth, for bond-moving (top) and the hier- 
archical lattice (bottom) for d = 3, divided by the initial cou- 
pling strength. The curves correspond to J = 0.3, 0.5, 1, 2, 4, 8 
from bottom to top, for both graphs. The data where ob- 
tained from a "cloud" of 50000 bonds. 

As mentioned in Section ||, the droplet picture predicts 
that beyond a length scale r ^ {J/hy^'^'^^^^^^ the contri- 
butions of the field and of the couplings to the free energy 
become comparable, and we expect that the strength of 
the couplings decreases beyond this scale. In order to 
test this prediction, we have plotted in Figure ^ the iter- 
ation number for which the two-spin coupling is largest 
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versus the logarithm of J /h. It should follow a law 



able y = qL^^'^ , we then obtain 



log2- 



ln2((d/2) -61) 



ln( J//i) + C ~ 1.15 ln( J//i) + C , 



with a suitable constant C. As the figure shows, the 
data for bond-moving agree nicely with this prediction, 
for the hierarchical lattice, the slope is larger and has a 
value around 1.4. This might be due to the fact that 
the field increases faster on the hierarchical lattice, lead- 
ing to an earlier reduction in the coupling strength than 
predicted by the scaling theory. 



.9 7 

2 fi 
5 
4 
3 
9 
8 

c 

■S 1 

2 R 
5 

4 

3 



























— \ ' 










^ ^ 


— \ ^ \ 


o - 



































2 3 4 5 6 7 

ln(J/h) 

FIG. 8. Iteration number for which the two-spin coupling 
J reaches its maximum, as function of ln{J/h), for different 
values of h. The symbols stand for h = 0.2 (triangle), h — 0.1 
(circle), h — 0.05 (square), and h = 0.02 (diamond). The 
lines have the slope 1.15 and offset C = 0.5 (bond-moving, 
top graph), and slope 1.4 and offset C — 0.3 (hierarchical 
lattice, bottom graph). 

Next, let us discuss the quantity A which is a mea- 
sure of the degree of non self-averaging. Figure || shows 
our results in the absence of a magnetic field. In the 
high-temperature phase as well as in the spin-glass phase 
A descreases with increasing system size and approaches 
zero, just as one can expect in the absence of RSB. At 
the critical coupling strength Jc, A remains constant with 
increasing system size, its value being A ~ 0.13 for the 
hierarchical lattice, and A ~ 0.15 for bondmoving. The 
constancy of A at the critical point can be explained from 
the scaling behaviour of the overlap distribution function 
P{q)- Critical scahng implies 



[Piq)] = L-^/^ \PiqL^/n 



and 



[P{q)P{q')] = L 



P{qLl^''')P{q'Ll^''') 



with (3 being the order parameter critical exponent, and 
v the correlation length exponent. Introducing the vari- 



(/y2F(y)dy 



independently of L. 
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FIG. 9. A as function of 1/J for /i = and L = 4, 8, 16, 32. 
The top graph is for bond-moving, the bottom graph on a 
hierarchical lattice. 

For low temperatures T = 1/ J , A seems to follow the 
law A ^ TL^^ with 6 ~ 0.24. This can be derived by the 
following argument: At low temperatures, most samples 
have a value of (g^) close to 1, and only a fraction p pro- 
portional to kTL^^ of all samples have system- wide exci- 
tations and have therefore some other value {q^) — x < 1. 
We therefore obtain 



1 — p -I- p[x\ 



and 



leading to 



[{qY]^l-p + p[x\ 



A'^p{l + [x^] - 2[x\) - kTL- 



In the presence of a magnetic field, we expect A to de- 
crease always to zero for large system sizes, because the 
system is always in the high-temperature phase without 
long-range correlations. However, as we will show in the 
following, A can become nevertheless very large for cer- 
tain combinations of the system size, the field, and the 
two-spin coupling strength. One can therefore easily get 
the impression that the system is not self-averaging, while 
in reality A increases only over a limited range of system 
sizes or parameters. 

Figure |l^ shows our results for A with a magnetic field 
h = 0.2. For J > Jc, the values are larger than without 
field, and they increase with increasing system size and 
decreasing temperature. We expect that as the system 
size increases further, A will reach a maximum and then 
decrease again. For fields stronger than h — 0.5, we see 
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this reversal in the trend of A already for the system sizes 
studied in the simulations. For J < Jc, Fig. shows that 
the curves for different L intersect each other, such that 
for high temperatures self-averaging is better for larger 
system sizes. Thus, the behaviour of A for weak fields 
seems to be qualitatively similar to that of the three-spin 
model. 

6 I ' , ' , ' , ' 1 
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FIG. 10. A as function of 1/J for h = 0.2 and 
L = 2,4,8,16,32 (from bottom to top for the 1/J = 0.5 
points). The top graph is for bond-moving, the bottom graph 
on a hierarchical lattice. 

For given system size and coupling strength J > Jc, 
there exists always a value of h for which A has a max- 
imum. The height of this maximum is higher for larger 
system sizes and for lower temperatures 1/J. Figure 
n^ shows the field for which A is largest as function of 
1/J. The data are in good agreement with a dependence 
hmax oc J, and hmax oc L~^'^^ for bond-moving. This 
means that A is largest when h '--^ JL^'^^'^^^\ For the 
hierarchical lattice, the fit to the data is best for a de- 
pendence h„iax C!C J, and hmax C!C L~°-^^. Just as in 
Figure ^, the effective value of d/2 — 6 appears to be 
larger on the hierarchical lattice than for bond-moving. 
We suspect that this is due to the fact that the field grows 
indefinitely on the hierarchical lattice. 




0.25 0.5 1 

1/J 



FIG. 11. The field value hmax for which A is largest, for 
L = 4, 8, 16, 32 (from top to bottom). The straight hnes are 
power laws hmax oc JL~^'^^ (top graph, bond-moving) and 
hmax OC JL~^'^^ (bottom graph, hierarchical lattice). 

These results can be understood if one considers the 
effect of the field on the overlap distribution P{q). With- 
out field, P{q) is a symmetric function, and varies con- 
siderably in shape for different samples for the system 
sizes typically used in simulations. This feature is seen 
in Monte-Carlo simulations ||2^ as well as in MKA [|. 
A magnetic field changes the shape of P{q) and moves 
the weight more and more towards positive q. In the 
limit /i ^ oo, all spins are aligned with the field, lead- 
ing to P{q) — 5{1). Along the boundary line L ~ 
(J/ft-)^/*^'^/^"^-', where the field is not yet strong enough 
to destroy all features of the low-temperature phase, we 
can expect that at least some samples still have large 
droplets than can be flipped without much free energy 
cost. In Monte-Carlo simulations ||lj], one finds indeed 
for certain intermediate parameter values a [/'(?)] that 
has a pronounced peak at some large q value, and a long 
and thin tail that extends almost all way down to g = — 1 . 
The authors point out that this feature results from most 
samples having only the main peak, and other samples 
having an additional second peak for some other value of 
q. They go on to argue that this is a non self-averaging 
feature characteristic of RSB, and that it should not be 
expected if there was no RSB. However, they also ad- 
mit that their simulations do not show a second peak 
in [P((7)] at a value qmim which is expected from Mean- 
Field Theory. Although we have not determined P{q) 
in the presence of a magnetic field within MKA, we can 
conclude from the behaviour of A that P{q) must have in 
MKA approximation exactly the same features that we 
just described for the Monte-Carlo simulations. Indeed, 
it is easy to show that A becomes large if most samples 
have a P{q) with one narrow peak at go, and some sam- 
ples have two peaks. For those samples with one peak, 
we have a small variance of q, 

Xs = {q^ -ql), 

which is essentially sample- independent. For those sam- 
ples with two peaks, we have a large variance xi which is 
different for each sample. If the fraction of samples with 
two peaks is p, we obtain 

^ ^ Ji^pki+MxfL „ 1 

{{I - p)Xs + p[xi]? 

p[x'i] 

where we have only kept the leading terms. As long as p 
is not much smaller than [x/]/Xsj ^ is of the order \/p. 
Thus, if Xs is small and p is small but not too small, A is 
large. The second condition is satisfied if the field is such 
that a small fraction of samples have a second peak in 
P{q), the first condition is better satisfied for larger L or 
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smaller T . This explains why we observe the maximum 
of A for those h values where the contribution of the field 
to the free energy is comparable to that of the couplings, 
and why the maximum of A is larger for larger systems 
and lower temperatures. Of course, for some even larger 
value of L, we expect P{q) to start having less sample-to- 
sample fluctuations, because the samples should become 
self-averaging. Then the argument will break down, and 
A should remain small. However, this range of system 
sizes seems to be beyond the reach of our simulations. 

In conclusion, we have shown that there exists a wide 
range of parameters over which the degree of non-self- 
averaging appears large for system sizes typically used in 
computer simulations. We expect our results to be valid 
beyond MKA. As we have argued for systems without 
a field , the apparent non-self-averaging must be at- 
tributed to the influence of the zero-field critical point. 
This influence reaches surprisingly far and creates a line 
in the h — J plane along which non-self-averaging is par- 
ticularly large, and which is somewhat reminiscent of the 
de Almeida- Thouless line. We are here in agreement with 
Huse and Fisher who argued already almost 10 years 
ago that Monte-Carlo simulation data for a spin glass 
in a magnetic field are strongly affected by the critical 
point. 



VI. DISCUSSION 



We have shown that for the three-spin model as well 
as for the spin glass in a magnetic field, a large degree 
of non-self-averaging found in computer simulations does 
not represent unequivocal evidence for RSB, but can be 
caused by finite-size effects. 

It seems, however, that a study of the non-self- 
averaging parameter A using Monte-Carlo simulations 
might be able to discriminate between RSB and the 
droplet picture for three- and four-dimensional spin 
glasses. As we have shown, A has in zero field a max- 
imum at Tc, and decreases again with decreasing tem- 
perature in MKA. If there was a low-temperature phase 
with RSB, the low- temperature value for A should prob- 
ably be larger than the critical value, and A should there- 
fore increase with decreasing temperature. Also, for tem- 
peratures below Tc, we found that A has its maximum 
not at zero field, but at some finite field value. The 
degree of self-averaging decreases deep in the supposed 
low-temperature phase. We expect a similar behaviour 
from the Monte-Carlo simulations. This would be a 
hint that non-self-averaging is strongest along the bound- 
ary between the field-dominated and coupling-dominated 
regime, and not in the region where one would expect a 
low-temperature phase with RSB. 
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